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HUMAN DYKi^MIC ORIENTATION MODEL APPLIED TO MOTION SIMULATION 


Work on this project has led to a Master’s Thesis by 
Joshua D. Borah. The thesis is abstracted below. 


Human Dynamxc Orientation Model Applied to Motxon. Simulatxon 


Joshua D. Borah 


ABSTRACT 


The Ormsby model of dynamic orientation, in the form of 
a discrete time computer program, has been used to predict 
non-visually induced sensations during an idealized coordi- 
nated aircraft turn. It was found that attitude and angular 
rate perceptions may be contradictory and furthermore, in a 
three rotational degree of freedom simulator, it is impossible 
to duplicate both simultaneously. To predict simulation fidelity, 
a simple scheme was devised using the Ormsby model to" assign 
penalties for incorrect attitude and angular rate perceptions. 
With this scheme, it was determined that a three rotational 
degree of freedom simulation should probably remain faithful 
to the attitude perception even at the expense of incorrect 
angular rate sensations. Implementing this strategy, a simu- 
lation profile for the idealized turn was designed for a Link 
GAT-1 trainer. Use of a simple optokinetic display was proposed 
as an attempt to improve the fidelity of roll rate sensations. 


Two open loop subjective tasks were designed, to obtain 
attitude and roll rate perception indications. A series of 
experiments were performed in our modified Link trainer to 
test the effectiveness of the tasks and to check model predic- 
tions and visual display effects. 



PROGRESS REPORT 


INTEGRATION OF ’VISUAL AND MOTION REQUIREMENTS 
FOR FLIGHT SIMULATION AND RIDE QUALITY INVESTIGATION 


June 197& through December 1976 


The followxng report briefly summarizes the work carried 
on during this period. Included are sections on visual cues 
dn landing, comparison of linear and non-linear ^?ashout filters 
using a model of the vestibular system, and visual vestibular 
interactions (yaw axis). One of the major accomplishments .of 
this period xias the completion of Hr. Joshua Borah’s master’s 
thesis, a copy of which is being sent separately. 




The subjective responses were self consistent, and. 
both tasks are considered to be useful for obtaining low 
frequency information. An unexpected difference was found 
between subjective indications and model predictions for 
the turn simulation. It can probably be explained by the 
response lag inherent in the task (low bandwidth) plus con- 
sideration of dynamic detection threshold effects; but this 
must be clarified by further work. The optokinetic display 
was found to be insufficient to significantly improve roll 
rate perception fidelity in the turn simulation, probably 
due to the short duration of the movements involved. 


Although not designed for the purpose, the predetermined 
simulation profiles were rated for realism by two pilots. The 
results did not contradict model predictions, although support 
was weak. A dynamic simulator motion logic was proposed, in- 
corporating the strategy derived from the model. Its use 
would enable pilots to "fly" the simulator, and may provide 
more convincing data for use in evaluating and revising the 
fidelity prediction scheme. 
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Work on this project has led to a Master’s Thesis by 
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Human Dynamic Orientation Model Applied to Motion Simulation 


Joshua D. Borah 


ABSTRACT 
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The Ormsby model\^f^ dynamic orientation, in the form of 
a discrete time computer program, has been used to predict 
non-visually induced sensations during an idealized coordi- 
nated aircraft turn. ItVas found that attitude and angular 
rate perceptions may be^' contradictory and furthermore, in a 
three rotational degree oft freedom simulator, it is impossible 
to duplicate both simultaneously. To predict simulation fidelity, 
a simple scheme was devisedi^using the Ormsby model to assign 
penalties for incorrect attitude and angular rate perceptions. 
With this scheme, it^was determined that a three rotational 
degree of freedom simulation should probably remain faithful 
to the attitude perception even at the expense of incorrect 
angular rate sensations. Implementing this strategy, a simu- 
lation profile for^the idealized turn was designed for a Link 
GAT-1 trainer. Us^e of a simple\^optokinetic display was proposed 
as an attempt to improve the fidelity of roll rate sensations. 

\ 

V 

Two open loop subjective tasks were designed, to obtain 
attitude and roll rate perception indications. A series of 
experiments were performed in our modified Link trainer to 
test the effectiveness of the tasks and to check model predic- 
tions and visual display effects. 



VISUAL CUES IN LANDING 


Introduction 

The overall aim of the research xs the development of 
practical tools ^/hich can extend the state of the art of moving 
base flight simulation for research and traxnxng. The immediate 
goal is the determination of the relative importance importance 
of various visual cues in flight simulation. The experiments 
to be conducted for this research are intended to obtain the 
perceptual response of humans to deviations from a nominal flight 
path during landing approaches. An existing fixed-based aircraft 
simulator is being modified to use a television projector with 
which the subjects will be shown recorded television Images of 
landing approaches. Verbal estimates of the magnitude of flight 
path deviations will be made hy the subjects, and this data will 
be used to construct a statistical model of the subject responses. 

Five subjects 'cjere run through the experiment to help 
refine the experimental protocol, evaluate the video tape con- 
figuration and quality, and provide preliminary data. The results 
are given later in this report. 



Physical modifications of the simulator room were 
completed. These include the Installation of the window 
shade and door shades, construction and mounting of the 
reflecting mirror for the projector, final positioning 
of the cockpit and minor rewiring of the electronics rack 
for a clearer field of view, and removal of the extra 
pane of glass from the cockpit window. 

Sideways shaking of the image was an annoying problem 
in the preliminary tests. It was at first thought to be 
due to mechanical problems in the Redifon at Langley, but 
was eventually discovered to be due to electronic incompat- 
ibility between the Amphicon. projector and the 1/2 inch 
video tape format. (The Amphicon has a 'slow' horizontal 
AFC instead of the more modern 'fast' AFC.) The problem 
was traced to the feedback loop around the horizontal oscil- 
lator in the Amphicon sweep chassis, and has been largely 
alleviated by modifying this circuit. The Amphicon will 
be realigned before running any more experiments . Further 
modification can then be done, if this proves necessary. 

Current configuration of the video tapes is as 


follows : 



(1) Long orientatxon run (backt-rards) 

(2) Long orientation run. nominal approach 

(3) Four scaling runs: minimum and maximum 

deviations at each 
distance 

(4) Long orientation run 

(5) Eight practice runs 

(6) Long orientation run 

(7) One hundred and tvrenty data runs 

Separate scaling and practice runs are provided for 
glide path and aim point tests. Criteria for the practice 
runs were that they have equal numbers of positive and 
negative deviations, large and small deviations, and short 
and long distances (but not all possible combinations) 
Also, deviations in the stimulus not being overtly tested 
were kept to the minimum level. All practice runs were 
presented in a random order. The original Langley tapes 
were edited and copied to make a set of tapes in the cur- 
rent configuration. 







The original tapes were recorded at Langley on a Pana- 
sonic NV-3020 video tape recorder. This type of recorder must 
be started and stopped by hand and has no electronic siritching 
to avoid recording start-up transients, so the control track is 
destroyed bet^reen runs. During playback, there is insufficient 
time for the player to re-stabilize itself after the loss of 
the control track, and the Amphicon projector cannot lock on 
to the signal before the start of the next run. The resulting 
degradation of picture quality, although lasting for only a 
second or so, is unacceptable for such short and carefully timed 
runs. Conventional video editing (such as is used for creating 
the scaling and practice runs) often makes the problem worse. 

In principle, "insert" editing should be able to solve 
the problem by recording the old video signal, independently 
of the control track, onto a tape with a brand-ne^j control track. 
This method assumes that the editor's playing deck can provide 
a usable video signal despite the damaged control track on the 
original tape. A reliably functioning video editor capable of 
performing good insei'ts (which require such refinements as 
flying erase heads) has not been available to date, so the 
insert method has not been properly tested. If this approach 
fails, new tapes will have to be made. New tapes of different 
landing approach conditions will eventually be needed any\jay, 
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and discussions of better methods of preparing them are in 
progress. 

When the original tapes were made, the skyplate on the 
Redifon simulator was not properly calibrated. The skyplate 
covers the lens to show a blank ’sky’, and is used to establish 
the visible ceiling by being partially opened. Lack of cali~ 
bration caused this ceiling to vary from run to run, but only 
one subject noticed this. It is not considered very important 
as long as the entire runway is visible, which is the case in 
all of the runs on the current tapes. 


Instructions to Subjects 

The purpose of this experiment is to determine your 
ability to detect errors in glidepath and aim point during 
aircraft landing approaches The experiment has two sets 
of video-taped landing approach runs (To save time, only 
a short segment of each run is shown.) During each set, 
you will be asked to estimate either glide path or aim 
point errors for each run. Both kinds of errors may occur 
simultaneously, but you should estimate only the one asked 
for. Since altitude along the glide path and aim point 

/ 
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miss distance depend on initial distance from the runway, you 
should base .your estimates on the ANGLES ojf the glide path 
and aim point vector errors (see description below) . 

Each set of runs begins ’5?ith two orientation runs to 
show you the toudhdown point and a correct approach. Four 
scaling runs follow to show you the largest errors in that 
set for either glide path or aim point. You should call the 
maximtim positive and negative errors "+10" and "-10" respec- 
tively, and estimate all other errors in terms of the -10 to 
+10 scale. For example, a positive error half as large as 
the maximum should be called +5- Except for orientation runs, 
there are no normal approaches (with error equal to zero) . 

The "glide path" is the path through space that would 
take you to the runway touchdom point. The correct glide 
path is the "glideslope" , which here makes a 3° angle to the 
horizontal. 

For any given glide path error, the difference in 
altitude will change with the distance from the runway, so 
you should estimate the ANGULAB error of the glide path (the 
glide path error angle) . 
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The "flight vector" is the direction you are moving 
in through space. The "aim point" is the place on the 
ground that you will reach if you continue along your present 
flight vector. The correct aim point is simply the runway 
touchdown point; to reach it, the flight vector must exactly 
align with the glide path. 

In an actual aircraft, only the instantaneous flight 
vector angle can be controlled directly, not the ultimate 
aim point, and this experiment is set up accordingly. For 
any given flight vector angle error, the ultimate touchdown 
point depends on the initial distance from the runway. Also 
the absolute size of the aim point error is not S3m!metrical 
for initial symmetrical flight vector angle errors. So you 
should estimate the error of the flight vector AEGLE^ rather 
than the ground distance of the resulting aim point. 

Note that it is possible to reach the correct touchdown 
point, even if the glide path is incorrect, and that the aim 
point can be in error even if you start out on the proper gJide- 
slope. If the flight vector is not aligned with the glide path, 
you may notice a slight change in the glide path during the run. 
If so, simply estimate the average glide path (or aim point), 
or that at the middle of the run. 

// 



A score of your performance during the test will he kept. 
You x^lll not be scored on correctly estimating the exact size 
of the error, just the right direction (+ or . Your score 
is simply the total number of estimates in the right direction. 
Your score does not represent your actual abilities as a pilot 
in a real aircraft and X'7ill be kept confidential. 

The runs average about 8 seconds long each, with 3 seconds 
between runs, so you should make your estimates quickly. You 
will have 8 practice runs, and you may repeat the scaling and 
orientation runs if you wish. 

Method and Preliminary Results 

The current experimental setup is as follows : 

Glidepath; 3° nominal 

±0.5°, ±1.0° deviations 
(±1.5° training) 

Flightpath* 3° nominal 

±0.6°, ±1.2° deviations 
(±1.8 training) 

3000 ft, 6000 ft 

// 


Distance: 



The glide path xs the vector from the aircraft to the 
desired touchdown poxnt on the runway. The flightpath is 
the velocity vector of the aircraft. 

Subjects are asked to give verbal estimates of the 
magnitudes of deviations on a subjective scale of -10 to +10, 
corresponding to the maximum deviations seen during the train- 
ing runs. Each subject estimates either glide path or flight 
path, but not both, during any particular experiment. 

A statistical model of the responses is constructed along 
the following lines: 

Response = Mean + (Glide path) H- (Flight path) 

+ (Distance) + [...interactions...] 

+ error 

The model is not intended to establish any cause and effect 
relationships, but to establish instead the relative importance 
of the different visual stimuli. 

For each experiment, we have the following stimuli: 

4 glideslopes x 4 flightpaths x 2 distances 
X 3 replications =96 stimuli 


Each subject sxts through two experxments, one for glxdeslope 
and one for flxghtpath. 


Plots of some of the results from a preliminary set of 
experiments are shown in the following figures. The actual 
magnitude of the deviation stimulus is given by the horizontal 
scale; the subject’s response to that stimulus (as derived 
by the statistical model) is given by the vertical scale. 

The sigma bars, where present, represent ±1 standard deviation 
of the "error" in the response which could not be statistically 
explained by the model. 

Keep in mind that the plots represent partial data from 
preliminary experiments. The subjects had different backgrounds 
of flight and simulator experience and variations in the experi- 
mental protocol xrere tried out. The plots were used primarily 
to test the application of the model, and do not all correspond 
to the same statistical confidence level or include all signi- 
ficant interactions. Nevertheless, they should give some indi- 
cation of what the final data may look like. 
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COMPARISON OF LINEAR AITO KONLINEAR WASHOUT FILTERS USING 


A MODEL OF THE HUMAN VESTIBULAR SYSTEM 


The purpose of thxs portion of the research is to 
discover Tjhat information the model of the human vestibular 
system can give about the effects of washout filters on a 
pilot’s perception of motion. A qualitative evalutation of 
linear versus non-linear washout filters as used in an air- 
craft simulator shows that the non-linear filter provides a 
"better" representation of actual motion than the linear 
filter. This subjective pilot response seems to be due pri- 
marily to the fact that the non-linear filter eliminates the 
false rotational rate cues presented by the linear filter. 
Hopefully, the model of the vestibular system will allow 
more objective results to replace the rating of "better", 
thereby providing a more quantitative method for predicting 
simulator performance. 

The model of the vestibular system, in its present form 
as developed in the Man-Vehicle Laboratory, exists as a FORTRAN 
program i^ith a very limited input/output format. The data, as 

AO 





recorded during the subjectxve tests at Langley, does not con- 
form to that format, so the following changes are being made 
to that program: 

(1) The program must, for convenience, run on the 
Laboratory’s PDP 11/34, Since it was prevxously 
run on the IBM 370, thxs requxres only a few 
modificatxons to the code. These changes were 
made and the program xs no\7 operational on the 
PDP 11/34. 

(2) The xnput format must be made more flexible, 
not only for the sake of thxs project, but for 
future projects as well. The program is in- 
flexible mainly because the state transxtxon 
matrices and Kalman gaxn matrices are pre- 
calculated for a state transitxon matrxx update 
interval of 0.1 second and a Kalman fxlter update 
interval of 1 second. If different update inter- 
vals are desired, these matrix values must be 
recalculated before the program is used. 

This continued recalculation doesn’t make much 
sense iii light of the new uses of the program. 
Therefore, the program will be changed to accept 
as input the two update intervals . Then the 
program t/ould calculate the matrix values and 



proceed as before. In order to accomplish this. 
It IS necessary to determine how these matrix 
values were first calculated. Then a subroutine 
must be added to the program to calculate the 
values, given the update intervals. This has 
been completed. 

(3) The output, in its present form, consists of 

pages and pages of data describing the various 
responses of the vestibular system. Since the 
PDF 11/34 IS equipped with a wide variety of 
graphics capabilities, it is desirable to display 
this data on graphs, rather than by hand plotting 
of the printed data. This part of the conversion 
IS not yet complete. 

Once the graphics changes are made to the program, it 
will be ready to interface with the data, and the analysis of 
the results can begin. 




VISUAL-VESTIBULAR INTERACTIONS (YAW AXIS) 


Effort durxng this reporting period has concentrated on modelling 
rotational motion sensation dependence on 'visual and vestibular cues, 
with exclusive concentration on yaw rotation about an earth-vertical 
axis. The primary results discussed below include, a) the proposal 
of a velocity drift model which helps explain low frequency drift 
observed during certain visual field presentations; b) a parallel 
channel visual- vestibular model which extends the earlier results 
obtained from combined input experiments; c) a human operator model 
appropriate to the velocity nulling task used to measure sensations; 
and d) a remnant correction technique useful for obtaining a more 
accurate estimate of the human operator’s linear characteristics. 

This progress report is divided into 7 sections . Sections 1 and 
2 discuss velocity drift results obtained in a velocity nulling experi- 
ment obtained under two different visual field presentations: a field 

fixed with respect to the observer, and one moving at constant velocity 
V7ith respect to the observer. Section 3 then presents some very simple 
functional models consistent \Jith this observed drift behavior. To 
extend the visual-vestibular model reported on in the May progress 
report, section 4 proposes aiid describes a dual-input parallel channel 
model, and develops the equations necessary to extract the model charac- 
teristics from the experimental data. Section 5 describes some of the 
results obtained from a dual-Jnput velocity nulling experiment, and 
presents Bode plots defining the human operator's combined controller/ 
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estimator describxng functions, appropriate to a parallel channel 
structure. Section 6 then proposes how the control portion of the 
operator’s transfer function can be divided out of the dual- input 
results, and section 7 presents various estimates of the inferred 
estimator describing functions resulting from such an approach. 




1.0 VELOCITY DRIFT WITH A FIXED VISUAL FIELD 


Review of previous data obtaxned from the first two series of 
closed-loop velocity nulling experiments and the results from a third 
series of similar experiments suggest that the velocity drift charac- 
teristics reported in the May progress report be amended. Specifically, 

It should be noted that: 

a. The latencn es and dual ramp drift characteristics illustrated 
in figure 2C of the May progress report should be discounted, since they 
may simply be an artifact of the drift measurement technique applxed to 
the strip chart histones. It has been found that computed values are too 
sensitive to slight variations in fitting the data in a piecewise linear 
manner. A more conservative approach has since been taken, in T’liich only 
one straight line is used to approximate the drift rate resulting from 
each fixed field presentation Such an approach is considerably less 
sensitive to variations in the fitting procedure, and has been reapplied 
to all of the presentations obtained in the second series of experiments. 

b. The drift rate statistics given in figure 2C of the May progress 
report should also be discounted, since they were obtained by pooling the 
results of a relatively well-controlled experimental series ( Series II) 
with those of a pilot series (Series I ), and the latter did not incorporate 
an adequate balance for presentation order, fatigue, and learning the task. 

The results presented here were obtained from two experimental series. 

^ One has already been described in the May report (series II) and was the 
source of the describing function (DF) data presented there. It suffices 





to note here that the subject used a spring-centered stxck for sub 3 ect 
control of velocity. The other series (series III) was a duplicate, 
except that a control wheel afforclxng no centering cues was substituted 
for the stick. Thus, the effect of possible centering cues could be 
investigated by comparing the results of the two experimental series. 

Both experiments were properly balanced as to the order of stimulus 
presentation, and the drift measurements were obtained by using a single 
straight line to fit the data. 

Both series used the same six subjects: four received 2 FIX 

presentations (that is, a presentation of a visual field fixed with 
respect to the subject, lasting 128 seconds), and two subjects received 
3 FIX presentations. Thus, in each series, there were 14 opportunities 
to observe possible velocity drift in the subject responses. 

In the first series using the control stick, there were 8 instances 
of observed drift, and 6 no responses (NR) observed, within the accuracy 
afforded by the strip chart recording. Counting each NR as a 0.0^/s^ 
drift rate, the population statistics are given by. 

Stick control: = 0.004°/s^j « 0.041°/s^; = 14 (la) 

A t-test shows that the mean is not significantly different from zero 
(p > 0.5), which IS what would be hoped for, since a significant non-zero 
mean would suggest a directional bias in the experimental equipment, pro- 
cedure, and/or subject population. 

We can now ask the following question. Is the self-centering property 
of the stick affecting the measured population response by providing a cue 
as to where the null is? Or, stated differently, are the subjects using 
the stick centering cue to augment their perception of the low-frequency 


motion cues’ 





With the experimental series repeated using the non-centering wheel 
as a control element (senes III) , there were again 14 opportunities to 
observe velocity drift. Actually observed were 13 cases of drift and 
one no response (NR). Again, counting the NR as a 0.0°/s^ drift rate, 
the population. statistics are given by: 

V7heel control: U2 " 0.015°/s^; = 0.050°/s^; = 14 (lb) 


As with stick control, a t-test shows that the mean drift rate is not 
significantly different from zero, suggesting the absence of a directional 
bias. Of more interest, however, is the question concerning a different 
population response due to wheel versus stick control. Comparing the 
statistics of (la) with those of (lb), we find that an F-test on the 
variances shows them not to be significantly different (p > 0.2), so that 
we can pool them (o^ = 0.046°/s") and use a t-test on the means. We find 
that the means are also not significantly different (p > 0.5), so that 
this statistical measure shows no difference between wheel and stick control 
Perhaps, however, this conclusion is biased by the fact that the NR' s 
of each series were included to arrive at the means and variances of (1) . 

By excluding them, and simply looking at the drift statistics of the 
responding population, we find the following: 


stick control: = 0.016; cf^ = 0.052; ^3^ ~ 8 (2a) 


wheel coirtrol. ^2 ~ 0.006; CT2 “ 0.056; N2 = 13 


(2b) 


An F-test shows the variances to be not significantly different (p > 0.2), 
so that we can pool them (o^ = 0.054 /s ) and use a t-test on the means. 
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Again, we find thatthe means are also not significantly different (p > 0.5) 
so that even excluding the NR’ s from the data, we find no significant dif- 
ference between vjheel and stick control, by these measures. 

Ilhat should be obvious at this point, however, is that the number 
of NR’ s observed with wheel control (1) is quite a bit smaller than the 
number observed with stick control (6) . To test the significance of 

this observation, we use a contingency table and the X^~test; 

0 

Stick Wheel 



Thus, there ^ a significant difference between stick and wheel control, 
in terms of the number of times zero drift (NR) was observed. The sug- 
gestion is that the stick provides centering information which completely 
suppresses drift in some cases, although the average drift rate is in- 
dependent of the type of control used. 

Since the statistical tests done above on the means and variances 
of tbe drift rates showed no significant differences between stick and 
wheel Control, it seems reasonable to pool the data. Of interest, then, 
is the manner in which tbe NR's are handled. If we assume the one NR 
observed with wheel control is a legitimate case of zero drift, uncor- 
rupted by a controller centering cue, then we are obliged to include it 
in the population results. This Is not unreasonable since it seems 
safe to assume that no controller centering cues were possible with 


the ^7hcel control. 



Turning now to the KR’s observed with stick control, one approach 
IS to simply exclude them all, on the basis of possible response cor- 
ruption due to centering cues. The corresponding contingency table 

test results in a value of 0.39, a considerable reduction from the 

0 

A. 76 value obtained above, and suggests that this is the proper direction 

in ^ifkich to proceed. Including only one of the NR’s observed during stick 

control results in a value of 0.11, and including too NR’s results in 

0 



value, thus, the minimum value is obtained T7ith one NR included in 

0 

the stick responses. 

For the results of the two experimental senes to be most congruent, 

in terms of NR occurances, it is clear that a contingency table test should 

result in a minimum value. Thus, the decision may be made to eliminate 

0 

5 of the 6 NR*s obtained under stick control. When the data is so edited, 
keeping one NR from each series, the following statistics -result : 

Drift rate; p = 0.011®/s^; cr = 0 050‘^/s^, N = 23 (3) 

A t-Lest shotjs the mean to be not significantly different from zero (p > 0.5), 
as expected. 

To gam an appreciation for the magnitude of the drift rates observed 
under fixed— field presentations, we can look at the statistics of the 
absolute values of the pooled sfcick/wheel data. 

Drift rate (magnitude): p = 0.043'^/s^; a = 0.027°/s^; N - 23 (A) 


These drift rate magnitudes are well below accepted threshold values for 
yaw axis earth-vertical rotation (- 0.30^/s^) and thus are consistent with the 




notion that the subject is completely unaware of his drift acceleration 
when deprived of vxsual motion cues. 

A summary of the above results is presented in Figure 1. 

A final note concerning the statistical characteristics of the 
velocity drift rates concerns the normality of the pooled data. Sho’s-m 
in Figure 2 is the cumulative frequency distribution (CFD) of the drift 
rates normalized with respect to the mean and variance of (4); super- 
imposed on this experimentally derived curve is the CFD for the unit 
normal distribution, N(0,1). Use of Kolmogorov-Smirnov test for normality 
strongly rejects non-normality, so that It is not unreasonable to conclude 
that velocity drift rates, caused by visual motion cue deprivation, are 
normally distributed. This will be contrasted to the results presented 
helo\r, concerning drift rate distribution observed with subject perforniance 
during presentation of a constant velocity (CV) visual field, 

A functional model of angular velocity perception incorporating 
the above-discussed drift characteristics will be presented in a later 
section, after discussion of the CV results. 




I ‘ VELOCITY DRIFT UITII >10 VISUAL CUDS 



© SIX SUBJECTS, 14 FIXED FIELD PKESENTATIOilS OVERALL, 
TVJO CONTROL METHODS 


CONTROL METHOD ] DRIFT OCCURMJCES ! 

■E (Vs^) 

O- (°/s2) 

STICK 8 

. 016 

.052 

WHEEL 13 

.006 

.056 


@ ONLY SIGNinCANT DIFEDRENCE BETWEEN METHODS IS IN wa*'5BER 
OP DRIFT OCCURAIJCEo. STICK PROVIDES A CENTERING CUE 
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2.0 VELOCITY DRIFT WITH A CONSTANT VELOCITY VISUAL FIELD 


The effect of a constant velocxty (CV) visual field presentation 
on subject performance in the closed- loop velocity-nulling task has 
already been qualitatively described in the May progress report. Briefly, 
it vas found that relatively large velocity drifts \^ere observed during 
such presentations, and that dri'ft always was in the direction of visual 
field motion. Given below is a more quantitative description of these 
drift responses. 

During the course of the previously-described velocity-nulling task, 
a constant velocity (4°/s) right-moving peripheral visual field xjas pre- 
sented to the subject. Velocity drift measurements wfere made on the 
ensuing subject responses, using the single straight line fit noted in 
the previous section. As in the case of fixed visual field presentations, 
rwo series v^ere run; one using self-centering stick control (series II) 

I 

and the other using wheel control (series III) . 

Both series used the same six subjects, four received 2 CV presen- 
tations (that IS, a presentation of a visual field moving at a constant 
velocity with respect to the subject), and two received 3 CV presentations. 
As with the FIX presentations, there were thus 14 opportunities per series 
to observe possible velocity drift in the subject responses. 

Out of a total possible 28 occurences of drift, there were observed 
27, vri.th one case of severe disorientation and subsequent inconsistent 
and task-unrelated response. This case has been eliminated from the data 
base whose statistics are given belo\j; 



stick 

control: 

= 0.261°/s^; 

= 0.141 Is ; 


(la) 

wheel 

control: 

= 0.328°/s^; 

©2 = 0.265°/s^; 

Nj- 13 

(lb) 


As in the preceeditig seel ion, we can ask xf stxck control results 
in sxgnif icantly different subject responses from those seen wxtb wheel 
control. It should be clear that they are equxvalent in terms of not 
suppressing drift responses, this in contrast to the large numbers of 
HR’s seeii with stick control and a fixed visual field presentation as 
noted in the preceeding section. Iftiat remains then is to compare the 
statistical measures above. 

Comparing the statistics of (la) with (lb) , an F-test on the variances 
sho^7S them to be significantly different (p < 0.05). A Welch t-test, how- 
ever, shovjs that the means are not significantly different (p > 0.5), and 
one is thus motivated to pool the data for the two control methods, in 

spite of the varaance differences. The pooled data for CV presentation 

« 

velocity drift rates are then characterized by the following statistics. 

9 2 

drift rate(CV): h = 0,293°/s^; cf = 0.209°/s ; N = 27 (2) 

It should be recalled that since all of the observed drifts were in the 
same direction (following the direction of the visual field motion) , 
these statistics also apply to drift rate magnitudes. It is also appro- 
priate to recall that these statistics apply to the single stimulus 
magnitude of a 4°/s right moving visual field. 

A simple t-test on the pooled statistics of (2) show the mean drift 
rate to be significantly different from zero (p < 0.005); this is to be 
contrasted to the approximately zero drift rate seen across the population 
in response to a fixed-field presentation. Further contrast between the 



3 / 


responses to the tx70 visual field environments is provided by comparing 
drift rate magnitudes. From (4) of the previous section, fixed-field 
drift was characterized by: 

Drift rate magnitude (FIX)* U *= 0.043 /s'; <? = 0.027°/s , N = 23 (3) 

An F-test on the variances of (2) and (3) above show a highly significant 
difference (p < 0.001), as does a Welch t— test on the means (p < 0.005) 
Thus we are led to conclude that drift rate magnitudes resulting from a 
CV visual field presentation are significantly different from those seen 
during a FIX presentation. 

Also of interest is the fact that the mean CV drift rate of approxi- 
o 2 

mately 0.3 /s is near three times the accepted yaw axis earth-vertical 
rotational acceleration threshold, suggesting a strong modulation of 
"vestibular" thresholds by visual inputs. A qualitative discussion of 
how such a drift rate can arise and remain undetected by the subject is 
given in the May progress report; a slightly more quantitative functional 
model is presented in the next section. 

A summary of the above findings is presented in Figure 1 

A final note concerning the statistical characteristics of the 
velocity drift rates concerns the normality of the pooled data. Sho\7n 
in Figure 2 is the cumulative frequency distribution (CFD) of the drift 
rates normalized with respect to the mean and variance of (2), presented 
in the same format used previously for the FIX drift rates. As before, 
use of the Kolmogorov-Smirnov test for normality requires us to reject 
non-normality (p > 0.2). Comparing this figure with the one drawn for 
the FIX drift rates (Figure 2, last section), however, suggests that the 
CV drift data is "less normal", because of the late rise and s3ow tail off 




I * VELOCITY DRIFT WITH CONqTiUiT VI^LOCIYY VISUAL FIELD 


© SIX SUBJECTS, 14 CV FIELD PRESEIJTATIONS OVERALL, 
rjO C015TR0L METHODS 


CONTROL METHOD 

DRIFT OCCURANCES 

^ (°/s2) 


STICK 

X 

13 

.261 

.141 

WHEEL 

14 

.328 

.265 


*one severe disorientation eliminated from talley 


© SIGNIFICANT DIFFERENCE BETWEEN TiTO METHODS IN VARIANCE, 
NOT ME2iNS 

G POOLED DATA- 

^ = 0.293 °/s^ O’ = 0.209 N= 27 


© MEAi^ SIGNIFICANTLY DIFFERENT FROM ZERO (P<.005), 

AJID ABOVE TillUiSIIOLD (a/ 0.1 °/s^) 

© RESULTS SIGNinCAtJTLY DirrElH:NT FROM FIX PRESENTATIONS 

© DATA BEST FIT WITH A LOG“NORi4AL DISTRIBUTION 
I4AGNITUDE ESTIMATE IMPLICATIONS ??? 






•1 



of the CFD. We are thus motivated to look at the log of the data, and 
the normalized CFD is plotted in Figure 3. Comparing this with that of 
Figure 2 shows that the CV drift rate is more accurately described as a 
log=normal random variable, rather than a normal random variable. How 
this IS to be interpreted is, at present, unclear, although it may sug- 
gest that some Weber la^7 estimation process is at ^7ork within each in- 
dividual's response, and is somehow reflected across the population 


response. 
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3 . 0 VELOCITY DRIFT MODELS 


The previous t\:o sections have discussed the statistical charac- 
teristics of velocity drift incurred by subjects during the task of 
closed— loop velocity nulling, when presented with two types of visual 
•motion cues. This section \?ill now present very simple functional 
models of the human operator which arc consistent with the observed 
drift behavior. 

3,1 Fi^ed Visual Field 

As discussed in the May progress report, past investigators have 
-attributed velocity drift arising out of a velocity nulling experiment 
(with visual motion cue deprivation) to be a characteristic primarily 
associated vrith the vestibular sensory channels. Functionally, this 
is an appealing interpretation because a very simple model can be 
constructed which is consistent with the observed drift and with the 
generally accepted properties of the vestibular sensors. Shovra in 
Figure 1 is such a model, a bilateral model of the two horizontal 
semicircular canals, whose outputs are differenced to provide an estimate 
of head angular velocity. Note that both canals are characterized as 
identical linear bandpass filters on velocity, but differing in DC gain 
and output bias. 

Shortly it will be sho^^n that a simple constant offset in the 
estimate w is sufficient to give rise to the drift behavior seen in the 
experiments. In particular. If a non-zero estimate (!) can he generated 
by the model of Figure 1, in the face of an angular velocity to which is 
actually zero, then it is a fairly direct matter to predict velocity drift 
in the closed-loop velocity nulling task. Of interest now is to see ho\? 





an offset can arise from a model parameter imbalattce". • From the figure, 
the angular velocity estimate is given by. 


w(s) = (bj^ - bj^) + 




w(s) 


( 1 ) 


(t^s-H) (t^s+I) 


It should be clear that if the bias terms are equal = b^^) , then a gain 
imbalance ^ not sufficient to provide a non-zero w when oj is 

actuallj'^ zero. Neither is it necessary, since a non-zero to can be simply 
obtained i;7hen w is zero by having the biases unequal (b^^ 4= b^^) . Thus a 
gain imbalance is irrelevant to a discussion of a steady offset in the 
velocity estimate, and for simiplicity mq. can assume a gain balance: 

V2 (2) 


so that (1) simplifies to the standard "cyclopean" canal model, with 
bias : 

K s 

oi(s) - to, (s) + w(s) I (3) 

(t^s+1) (t^s-I-I) 


where the bias velocity is defined by 


“b - V - \ 

The functional model corresponding to (3) and (A) is sketched in 
Figure 2, and v^'ill be used in the sequel. 

Now to sho\7 how such an angular velocity bias can give rise to 
the drift observed in the velocity nulling task, it is necessary to 
refer to the loop diagram previously presented in the Hay progress 
report and repeated in Figure 3. Since we are presently considering 
the purely vestibular situation, in wh^ch the subject is presented 
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with a stationary visual surround, the estimator block E(s) can be 
replaced by the vestibular model of Figure 2. Furthermore, as will be 
discussed at greater length in section 6, the control strategy 
C(s) can be approximated by an integrator with gain, or; 

C(s) = K /s (5) 

c 

Simple block diagram arithmetic then yields .the following expression 

for the angular velocity, w, resulting from the sub 3 ect*s remnant, n, 

his vestibular bias, w, , aiid the input disturbance, d: 

b 

= ThIcF 

where E’ is the linear portion of the vestibular estimator: 

E' (s) E K^s/{(t^s+1) (t^s+1)} (7) 

Since we are interested in low-frequency behavaor (specifically, drift), 

note that* 

lim E' (s) = K s (8a) 

s -^0 

Further, the trainer transfer function is unity at DC, so that 

lim G(s) = 1 (8b) 

s -^0 

Since thedisturbance d(t) is a sum of sinusoids (as described in the May 
progress report) , 

lim d(s) = 0 (8c) 

s-^0 

Finally, if we assume the remnant to liave no po\jer at zero frequency, then 


lim n(s) = 0 
s ->0 


(8d) 



Substituting (5) and (8) into (0), we fand the low frequency portion of 
the angular velocity signal to be given by: 


—fvlS. 

liin ci(s) = lim — )~T (9') 

s^O s-^0 ^^^c\ ^ 

where we have used the fact that the bias is a constant over time, 
so that 


to, (s) = 03, /s 
b b 


(Be) 


In the time domain, then, ( 9^ implies that, due to the velocity estimate 
bias tOj^, the subject will continue to accelerate at a constant rate, his 
angular velocity being given by; 


w(t) = - 


KK 

c 

1+KK Kt ‘^b*^ 
c 1 


(PIX Drift) 


( 10 ) 


assuming zero initial conditions. The minus sign, of course, implies 
that a positive (rightward) bias will give rise to a negative (leftward) 
drift. Thus, the simple cyclopean vestibular model, modified with the 
addition of a bias on the output, appears to be an adequate descriptor 
of subject performance when one is deprived of visual motion cues on the 
velocity nulling task. 

Presumably, the estimator bias for an individual could be inferred 
from (10) by measuring the drift acceleration and by estimating the indi- 
vidual transfer function gams K, K and K . Tlie same cou3d be done for 

c X 

the entire population, hut it is simpler to recognize from (3) of section 1 
that the average drift rate w is zero over the population. Thus, from 
(IQ) , the average velocity bias over the population must also be zero, 
since KK has always been found to be non-zero (see section 6). Estimation 

4 ^ 



of the variance of o), from drxft rate measurements is, of course, 

b 

complicated by the unkno^m statistical properties of the gains in 

( 10 ). 

One final qualitative aspect of this bias model is worth noting, 
and concerns the subject’s perceived velocity while engaged in the 
velocity nulling task. His estimated velocity to is neither zero nor 
(Oj^, but an intermediate value found by substituting the transform of 
(10) into the cyclopean canal model of (3), to yield; 


to K s -KK ^ 

T (t s+1)(t„s+1) ^l+KK K ^ I ^ 

s-^0 i 2 cl 

where the subscript LF indicates that we are interested in the low 

frequency portion of perceived velocity. The above expression slmplafies 

to yield the follo'oing relation between bias velocity and perceived 

velocity: 



1 

1+KK K, 
c 1 




( 11 ') 


Thus, the subject perceives that he is moving at a constant velocity, 

and hence is obligated to provide a compensatory control stick deflection 

to null it, thus leading to the eventual acceleration drift seen in the 

records. Note that his perceaved velocity co is simply the bias velocity 

LF 

attentuated by the closed loop gain, so that individuals with higher loop 
gains ^’ill tend to have lot^er values of perceived velocity, and vice versa, 
assuming other factors remain equal between subjects. 

3.2 Constant Velocity Visual Field 

To this point, we have been concerned ^jith the drift incurred with 
a fixed visual field (FIX); a similar argument can be made to heip under- 
stand the cause of drift under constant velocity visual field presentation (CV) . 





Shown in Figure A is perhaps the simplest possible parallel channel 
visual-vestibular model, in which it is presumed that visual surround 
velocity contributes only to the low frequency portion of the angular 
velocity estimate, in a linear manner. Some 3Ustification for this 
model has already been given in the May progress report; more will 
be given in a later section describing an experiment whose goal is 
to determine the components of such a parallel channel model. For now, 
it suffices to note that the velocity estimate from this model is 
given by: 


K^s 

fii(s) = to (s) + K-to„(s) + to (s) (12) 

^ (T^s+1) (T^s+l) 

■sphere (0^^ and to^ refer to vestibular and visual field velocities -res- 
pectively. Note that by defining 

+ (13) 

the equation is identical to (3), the biased estimator for the purely 

vestibular situation. Thus, the same low-frequency derivation is 

IS applicable and (10) can be used to describe the CV-induced drift 

rates seen experimentally, with (O^ in the equation replaced by (O^ above, or 

KK 






In the previous section it was noted that all CV- induced drifts 


were observed to be in the same direction as the stimulus visual field 
velocity. What this suggests is that the visual field effect is large 
with respect to the (bilateral) vestibular offset term. This is seen 
fairly directly by recasting (10) and (14) in terms of drift acceleration 
levels : 

(15) 

^CV ^ 

where A is defined by 




H-KK K, 
c 1 


(16) 


But, from the previous two sections, the average drift accelerations ■^yere 


^FIX 

0.04°/s 

"cv 

0.29'^/s 


(17) 


So that inspection of (15b) would lead one to conclude that, over the 
population, 

Wb « ^ 2^2 


that IS, the vestibular bias velocity is small with respect to the CV— 
induced velocity sensatxon. It should be recognized that this conclusion 
is applicable to the particular value of used in the experiment, a 
4^/s right-moving visual field, and lower field velocities may not allow 
similar conclusions to be made. 




4.0 PARALLEL CHANNEL VISUAL-VESTIBULAR MOUEL 




The previously described velocity-nulling expci*iments (series I, 
II, and III) looked at closed-loop velocity control in the face of a 
vestibular disturbance combined with one of three types of visual 
surround environments; a counterrotating visual field (CON) which 
exactly confirmed vestibular inputs, a fie]d stationary with respect 
to the subject (FIX) which provided no visual motion cues, and a 
field moving at constant velocity with respect to the subject (CV) 
which induced circularvection sensations. 

The results, analyzed in both the time and frequency domain, 
support the notion of a frequency separation of visual and vestibular 
inputs, where the visual cues provide the steady state or low frequency 
cues, i7hile the vestibular cues provide complementary high frequency 
or transient information. The human operator descrDbing function was 
calculated for each of the three visual field conditions and was 
modelled as a lag-lead function whose parameters were dependent on 
the particular field condition (see May progress report). 

The describing function itself relates actual trainer motion to 
the subject's control stick/wheel output, and thus can be viewed as a 
"vestibuiar" transfer function whose parameters depend on the parti- 
cular visua] motion cues being presented to the subject. This, of 
course, assumes that the subject's control strategy for nulling per- 
ceived motion is essentially unity throughout the frequency range of 
interest; this subject will be addressed later in section 6. 




In order to look more closely at \jhat Is essentially a dual~lnput 
problem, a parallel channel estimator model was proposed, and another 
experimental series (IV) initiated to see xf such a model could explain 
in greater detail how visual and vestibular inputs are combined to 
arrive at a sensation of motion. The approach is to work x;ith two 
describing functions: one relating trainer velocity to wheel output, 

and the other relating visual field motion to wheel output. 

In the previous three series, noise was injected into the trainer 
loop so that the subject would be obliged to provide compensatory 
stick or wheel deflections to null his perceived angular velocity. 

In the present dual- input (DI) senes, noise is also injected into 
the visual field loop, so that any visually- induced motion sensation 
must also be compensated for by the subject. By choosing the two 
disturbances to be independent of one another, it is impossible for the 
subject to perform both nulling tasks at once. By examining the results 
in the frequency domain, it then becomes possible to see what components 
of the subject's response are due to vestibular inputs, and what are due 
to visual inputs. From this, one can infer a DI describing function 
model of perceived motion. From experience with the previous series, in- 
jecting noise into the trainer loop (vestibular path) Is fairly straight 
forward: simply sum it ^Jith the subject's wheel signal prior to 

being sent to the trainer drive circuitry. The same approach can be 
taken for the projector loop (visual loop) : sum a second noise signal 
with the wheel signal and send this combined signal to the projector 
drive circuitry. 




A functional block diagram of the overall system, including a 
conjectural dual-input model of the human operator, is shown in 
Figure 1 . The two inj ected_ noise signals are denoted by d^^ and 

d^j and the wheel signal X is shown as an input to both the visual 

and vestibular loops For this series and others following, 
wheel control vzas used to avoid possible centering cues provided by 
the spring loaded stick. 'I\7o points should be noted First, the 
sign of the wheel signal is changed ^?hen it Is sent to the projector 
drive, to make the resulting visual field motion consistent \ 73 th the 
trainer motion. Thus, a right wheel motion results in a right trainer 

motion and a left visual field motion; i.e. analogous to the counter- 

rotating field situation of the previous three series. The second point 
to note is the addition of a prefilter, F^, m the visual field path, 
necessary because the projector drive alone, G^, has a relatively hagb 
bandpass compared to the bandpass of the filtered trainer transfer 
function, F2^1 choosing so that 

Vi ■ ''2S 

then, in response to the subject’s wheel deilections, the visual field 
motion will mimic the trainer motion, exhabiting the same ampJitude 
attenuation and phase lag over the frequency regime of interest Another 
way of saj’ing this, is, that, in the absence of any noise injected into 
either loop, the visual field motion, in response to the subject's wheel 
deflections, should be indistinguishable from the counterrotating visual 
field motion used in the previous experimental series. Thus, the pre- 
filter F 2 helps make the visual field a more compelling motion cue. 
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A linear model of the subject, interlaced with the experimental 
apparatus, is diagrammed in Figure 2. Here, it is assumed that the 
visual and vestibular channels v;ork in parallel, with their outputs 
summed to provide an overall estimate of self-motion. Tins is the 
same approach used in the velocity drift modelling of section 3, 
but is to be contrasted with the single channel "vestibular" model, 
augmented by visually- induced parameter variations, derived from the 
results of the previously described experimental series. As before, 
an internal model of zero perceived velocity is assumed to be a 
set-point (ca^ = 0) upon which acts a linear control logic C(s) to 
generate the appropriate hand motion which drives the control wheel. 
A remnant signal provides for a source of subject output response 
which is uncorrelated x^ith either of the disturbance inputs, and 



The remainder of this section will now be concerned with the 
manner in which the estimator transfer functions and E^ may be 
solved for ih terms of the three loop inputs, d^, d^, and n, and 
the three loop outputs, and X, 

From block diagram algebra, it is a direct matter to show that 
the three outputs are determined by the three onputs in the following 
manner : 

X(s) = |(Kn(s) + KPCE^d^(s) - KPCE 2 d 2 (s)) (a) ' 

Wl(s) = |(Kn(s) + (1 + KPCE 2 )d^(s) - KPCE 2 d 2 (s)) (b) (1) 

W 2 (s) = “ |(Kn(s) - KI?CE^d^(s) + (1 + KPCE^'ld 2 (s) ) (c) 



where 


A = 1 + KPC(E^ + E2> 


( 2 ) 
















and where the transfer function dependence on the Laplace- transform 
variable s has been omitted for clarity. We can now make use of 


auto- and cross-power spectral density functions to solve for and 
£ 2 * Correlating with X and we have, from (la) and (lb). 


^ - KPCE,^. , - KPCE„$. . ) 

Xd^ A ndj^ 1 2 d2d^ 


(a) 


(3) 




Since the remnant is defined to he uncorrelated with the input dis- 
turbances, then 


$ , = $ , =0 

ndj ndj. 


( 4 ) 


further, if the experimental design is such that and d 2 are un- 
coxre] ated, then 


$ 

a^d2 


= 4> 




= 0 


( 5 ) 


Thus^ dividing (3a) by (3b) and using (4) and (5) to simplify the 
result, we find that 


Xdj^ KCE^ 

1+KPCEp 

In a siBillar fashion, it may be sho’wn that 
$ 

"^Xd2 KCE 2 

^t02d2 i+IOPCE^ 


(6a) 


(6b) 


Since the 3eft-hand-sidcs of (6) are computable from the measured outputs 
of the experiment, we define / 



so that substxtuting into (6) allows for a solution for CE^ and CE^: 


^ ^ l+P^a^a^ 

_ 1 

2 K 9 

l+P^a^a^ 


(a) 


( 8 ) 


(b) 


Note that the control strategy C is embedded V7lth the estimators E- 
and E^ , as is to be expected, and cannot be separated from them in 
this type of experiment. Separation of control from estimation is 
the subject of section 7; the present discussion will be 
concerned with estimating the composite functions CE^ and CE^. 


4.1 Disturbance Inputs 

One major aspect of the experimental design concerns the choice 
of the two disturbance inputs d^ and Basically, they were chosen 

to meet the following requirements 

a. Their frequency content must span the frequency range of 
interest: -0.01 Hz to -1.00 Hz. 

b. The high frequency poi^er should gradually taper off, so as 

to avoid the "shelf” power spectrum used in earlier experiments. 

c. The amplitude must be reasonabJy sized in terms of human operator 
threshold and maximum limitations on the equipment. 

d. The two signals should be uncorrelated, to satisfy (5). 



Chosen were two "pseudo -random" signals, each composed of a 
sum of sxnusoids. Each sxnusoidal term xs an xntegral prxme multiple 
of a base frequency of 1/128 Hz (0.00781 Hz), so that the total pcrxod 
of each sxgnal is 128 seconds, a perxod whxch xs unrecognizable by 
the human operator in thxs task. Formally, the signals are defined by: 

13 13 

d = S a sin n w„t d_ = E b sin m o) t 

1 tX lO 10 

X 1=1 1=1 

where = 2 tt/T, T = 128 seconds, and where the n^ and m^ are alternating 
prime numbers so as to avoid harmonic ambiguities and assure a zero 
correDation betx<reen the signals. 

The amplitude spectra for the tT?o signals are given in Figure 3 
and arc given in amplitude ratio (AE.) form, referenced to a base low- 
frequency magnitude of 1.2*^/s. As can be seen from the plot, the frequencies 
of the two signals are interleaved, and both follow the AR curve associated 
with a double lag-lead transfer function, given by 


s^H- 2c w s + 

_± 1_J 1_ 

()^ s^+ 2C w s + 03^ 

2 2 2 


where 


(03^, (o^) = (0.475, 0 150) Hz 
(C^, C^) = (0.707, 0.707) 


The lag and lead break frequencies were chosen to give a gradual trans- 
ition beUJeen the large low frequency amplitudes and the small high 
frequency amplitudes; the 20 dB ratio beU?een the two extremes was 
chosen from past experience with disturbance Inputs into the trainer. 












To avoid rapid start-up transients when the two disturbances begin, 
the signs of the amplitudes are alternated so that* 

sign(a^^^) = - sign(a^) ; sign(b^^j^) =-sign(b^) 

where 

sign(a^) = 1 ; sign(b^) = 1 

\ — 

The appearance of the resulting signals in the time domain is as sho^-m 
in Figure 4 . 


4 . 2 Plant Dynamics 

The other ma^or aspect of the experimental design concerns the 
choice of the plant dynamics. For computational simplicity, it would 
be desirable to maintain unity gain and zero phase lag for the entire 
range of test frequencies, but, as noted in the May progress report, 
the trainer transfer function, G^(s), exhibits a strong resonance at 
about 1.5 Hz, due to the mechanicaJ properties of the load and drive 
system. To avoid this, a pref liter was added, F^(s), so as to ensure 
linear operation and reliable velocity feedback information over a 
lower frequency range. The combined plant, looks like a unity 

gam second-order system, with a break at 0.90 Hz and a damping ratio 
of 0.70. 

It xras noted earlier that this kind of response necessitates the 
use of a prefilter, F^, in the projector drive circuit, to ensure that 
F^Gj^ = ^ 2^2 prefilter was implemented on the analog computer 

and is a second-order filter X7ith a break at 0.92 Hz and a damping ratio 

<^o 


of 0.70. 














To ensure that the traxner and projector transfer functions, with, 
thexr associated prefilters, were close approximations of one another, 
input-output testing \'as performed on each. The wheel signal was held 
at zero while the disturbance was sent into the trainer drive and the 
disturbance d^ was sent into the projector drive. The time histones 
of the t;-7o, along with the resulting trainer and visual field motion, 
are shown in Figure 5 • Note the high frequency attenuation in both 
channels and note that the visual field velocity is the negative of 
the input command d^. 

By taking the Fourier transforms of these signals, the transfer 
functions of the trainer and projector drive can be computed, Showm 
in Figure 6 are the computed amplitude ratios (AR) and phase lags ((|)) 
for both the trainer and projector system; superaraposed on the data are 
the AR and curves associated with the second-order transfer function 
given in the figure. This data thus substantiates the experimental 
condition thaf 




FiG^ 




n 




w - 5.65 rad/sec (0.90 llz) 
n 


? - 0.70 
n 


and thus visual field motion v?ill mimic trainer motion, in response to 
both disturbance and control wheel inputs. Note that this knoT/ledge of 
the plant dynamics is a prerequisite for solving for the estimator functions 
and E^, as reference to (8) will show. 






4 . 3 Computatxonal Approach 


At thxs point, it is appropriate to provide a brief description of the 
computational approach used in processing the experimental data to arrive 
at describing function estimates. 

The analysis given above was based on the use of auto- and cross-po\/ei' 
Spectral density functions, but, from a practical viewpoint, it is computa- 
tionally simpler to work with the Fourier transforms of the measured signals. 
To see ho’tr this is accomplished, it is convenient to regard the v;heel response 
signal as being composed of three component signals, as follows. 

13 13 

X(t) = Z sin(n w t + tn -) + ^ + ^ ) + X (t) (10) 

i=l 

By reference to (9), the first component contains those frequencies 
associated vrith d^, the second those associated 'i7ith and the third, all 
other frequencies ^7hich are integral multiples of the base frequency 

CO 

X (t) = Z X sin(il w„t + 0 ) (11) 

r i=l ^ ^ i 

where the set of is the set of all integers excluding the sets of 
and {m^}. It should be recognized that "^jith this definition, (10) is an 
approximation to the actual wheel signal, since non- integral multiples of 
the base frequency have been excluded. However, such an approximation 
is consistent with the discrete frequency resolution which results from the 
digital Fourier transforms used in processing the data. 




Recognizing that A is periodic, with a period T = 211/(0^ = 128 seconds, 
then It IS a direct matter to solve for the cross-correlation function 
between the wheel signal and the "vestibular" disturbance: 



T/2 

(t) = A(t)d (t+T) = i / A(t)d (t+T)dt 
^ -T/2 ^ 


( 12 ) 


Substitution of (9), (10), and (11) in the above relation, followed by an 
application of the well-kno^m orthogonality properties of sinusoids, results 
in the following expression for * 
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1=1 


a A . cos(n.ti) T -({).) 
1 li 1 0 li 


(13a) 


so that only the frequencies of d^ appear in . A similar expression 

may be found for • Af we follow the same procedure as with the xrhcel 

signal, and separate the trainer signal into three components, we have. 


13 13 

u (t) = Z u sin(n w t + ) + Z to sin(ra to t + ) + w (t) ''j(l^) 

1 1 li 3 0 li . , 2i 1 0 2i Ir 

1=1 1=1 


It follows that 


, (T) 

Vi 


13 

Z a to cos(n to T - ij; ) 
1 li 1 0 lx 


(13b) 


• 60 , again, we have only the frequencies of d^ in ^ Kow, from the 

definition of ct^ given in (7), and the sinusoidal composition of the 
crosG-correlat3on functions of (13), it follows that: 




(x=l,13) 


jA(n u )j 

|a (n.w )( = - — = X- /(O- . = I 7- ^ 

( 15 ) 

3 c a^Cn^w^) = ^ X(ii^Wq) - ^ 01 ^( 11 ^ 03 ^) (x=l,i3) 

Or, more compactly, 

= X (h^Wq) / to^ (n^w^) (x=l,13) (16a) 

Naturally, the same results are applicable to the defxnition o£ a ^ gxven in 
(9); 


a2(m^03Q) “ X(m^ WQ)/c02(m^WQ) (i=l,13) (16b) 

Thus instead of bexng involved i?xth the computatxon of cross-correlation 
functxons and their transforms to po\jer spectral densities, the idcntifD cation 
problem becomes one of simply input-output transfer function computation 
via (16), to define and That is, Fourier transforms may be made 

directly on the measured signals X, o)^, and and the complex algebra of 
(16) can then be used to specify 01 ^ andcc^ at the discrete input frequencies 
of d^ and d^. 

It should be noted from (16) that since and are not defined at the same 
frequencies, then the computations for CE^ and CE^ indicated by (8) cannot 
be made. I-Jhat is required is an assumption of continuity, in the frequency 
domain, of the transfer functions introduced by the linear model. This 
then allows for one to interpolate across frequency to obtain the needed 
cc^ and values. That is, one can interpolate between the values 

■^7 



defined at the n (o. frequencies to obtain values for a, at the m frc- 
quencies, and conversely for a^. With a^-and both defined in this 
manner at all 26 input frequencies introduced by both d^ and d^, (8) may then 
be used to define the describing functions CE^ and CE^, for all 26 frequencies. 
This vas the approach used In the data processing which followed FFT processing 
of the recorded signals. 

This then completes the discussion of the parallel channel visual- 
vestibular model and the dual- input experimental and computational approach 
for estimating the model’s separate visual and vestibular transfer functions. 




5.0 


DUAL- INPUT EXPERIMENTAL RESULTS 



As in the previous velocity-nulling experiments, the subject’s task 
in this series was to maintain his sensation of self-velocity as 
close to zero as possible, by appropriate motions of the control 
wheel, and by inferring self-motion from the continued application 
of visual and vestibular cues. After a familiarization period ’tJith 
the procedure and equipment, the subject performed one continuous 
run of velocity nulling ^7hich lasted for approximately eight minutes. 
During this run, the vestibular disturbance (d^ referred to in the 
last section) was continuously inputted to the trainer drive. The 
visual environment, however, alternated between Wo modes: a counter- 

rotating" field (CON mode) which provided accurate confirmation of 
vestibular cues, as in the previous series, and an independently movang 
field (DI mode), which was the result of a visual field velocity dis- 
turbance referred to in the last section) inputted to the pro- 
jector drive. Each presentation mode lasted for 128 seconds, and 
alternated with the other, so that there were two presentations of each 
visual field condition to a subject 


Scries A CON, DI, CON, DI 
senes B. DI, CON, DI, CON 


Half the test population of 6 subjects received series A, while 
the other half received series B, so as to provide a balance for 
fatigue and learning effects when averages are taken across the popo" 


lation. 



Sho\TO xn figure 1 are time histories of a portion, of a subject's run, 
showing the vestibular disturbance the trainer and visual field 
velocities and and the subject* s compensatory response X. The 
first portion (CON mode)' illustrates good velocity^nulling performance 
when the subject is presented with a counterrotating visual field, 
and specifically shows his ability to null out low frequency disturbances, 
presumably because of the corroborating visual motion cue provided 
The second portion (DI mode) illustrates poorer performance, especially 
with regard to nulling out low-frequency drift in the trainer velocity. 
Presumably, his low frequency stick response is primarily dedicated 
to nulling out the visual field velocity disturbance as evidenced by 
the contrasting lack of drift seen in the field velocity history. 

5.1 Frequency Domain Results 

More definitive observations on subject response during DI presentation 
can be made in the frequency domain. Sho\m in figure 2 are two wheel 
deflection amplitude spectra plots, superimposed on one another, obtained 
from one individual by transforming the recorded wheel history via an 
FFT program. The sample rate used vas 8 Hz, and the sample length \<ras 
128 seconds, so that each FFT performed covered one entire DI presenta- 
tion , resulting in the two spectral sets &ho\ni 

The circles identify subject response at the frequencies contained in 
the "vestibu3ar" disturbance signal the squares identify response 
at the frequencies contained in the "visual" signal and the dots 
identify the remnant; i.e. response at frequencies contained in neither 
disturbance input. At the "vestibular" frequencies (circles) , the 
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response can. be seen to follow the shape of the disturbance spectrum 
presented in the last section (figure 3), indicating an appropriate 
compensatory response to trainer velocity. This, of course, assumes that 
the phase is appropriate, a subject to he discussed later. In contrast, 
response at the "visual" frequencies (squares) shows a sharp drop off at 
low frequency, indicative of little attention being paid to the high 
frequency visual inputs in the task of nulling perceived self-velocity . 

A smooth curve has been drawn through the remnant response (dots) so as 
to provide a simple approximation to the remnant power, and will serve 
as the basis, in section 8, for a discussion of remnant power correction 
to derived operator transfer functions. For noxr, however, it suffices to 
note the general trend of remnant dominance of subject response as frequency 
increases, especially at frequencies greater than approximate!}' 0.5 Hz. It 
IS also appropriate to note the difference in signal-to-noise power (S/N) 
ratios betweezi response at "vestibular" frequencies and at "visual" fre- 
quencies. To calculate the S/N ratios, we square the amplitudes of figure 

t 

2 to obtain the power spectra. Then, by linearly interpolating between 
responses at the "vestibular" frequencies (circles), and integration over 
the frequency range shown, the total "vestibular" power in the wheel 
response can be calculated. A similar calculation yields tbe total "visual" 
power, while integration under the smooth remnant curve yields a figure for 
remnant or noise power. Ratios of these 'figures then result in S/N ratios 
for both sensory channels. Since figure 2 illustrates response for two 
DI presentations, we obtain tbe following two sets of S/N ratios: 



7^ 



( S/N ratios xn dB ) 

Obviously major difference is the 10 to 12 dB difference between 
response at the two sets of frequencies, indicative of the greater 
importance of vestibular cues in this nulling task, especially at 
the higher frequencies. This aspect of the response will be discussed 
in greater detail when remnant corrections to the response are considered. 

It should be clear that subject response during the CON visual field 
presentations of this DI experiment can be similarly analysed in the 
frequency domain. Since many of the results are similar to those 
already presented in the May progress report (for earlier experimental 
series) , a discussion of the CON results will be deferred to a later 
section. For now, we will continue to describe the results obtained 
from DI presentations. 

5 . 2 Population Averages for CE ^ ^ and CE ^ 

Plots similar to that of figure 2 could also be dra™ to illustrate 
the amplitude spectra of trainer motion and visual field motion 
and qualitative conclusions could be dra\m regarding the appropriatness 
of the subject's compensation. It is more direct, however, to simply 
apply the computational techniques of the last section to this frequency 





data, and arrive at estimates of the two describing functions, 
and CE„. This has been done for both DI presentations to each subject, 
for all SIX subjects, and the resulting population average Bode plots 
are given in figures 3 and 4. 

The data points in the figures identify averages for the six subjects, 
while the smooth lines drawn through them simply indicate trends as 
the frequency changes. Several points are worth noting First, the 
gain for the "vestibular" describing function, CE^, follo^JS, in the mid- 
range frequencies, what would be expected from a lag-lead function. It 
may be recalled from the May progress report that a lag-lead function 
formed the basis for the adjustable parameter vestibular model, and the 
results sho\m here support that approach. The earlier results also show 
a lead at high frequencies, again evident in figure 3. The major dif- 
ference, ho’i^ever, is in low frequency behavior* the earlier single channel 
model indicated a levelling off to some fixed DC gain, \diereas the gam for 
CE^, IS seen to be dropping off as the frequency approaches zero. The 
"washout" characteristic is entirely consistent with our notion of neg- 
ligible canal response at very low frequencies, because of the canal's 
AC physical properties. The Bode plots defining the "visual" transfer 
function, CE^s sIiovt quite contras Ling behavior. At low frequencies, the 
gam IS higher than in the vestibular channel, supporting the previous 
statements concerning the importance of lo\;-f requency visual cues in 
determining motion sensation. Up to approximately 0.1 Kz, CE^ looks like 
a simple integrator (in gam and phase) , vzhich, as ^7ill be seen in the 


next section, is simply a reflection of of the control strategy used by the 
subjects. Above that frequency , the CE„ gam levels out, followed by a slight 

t ^ t 

















lead at the highest frequencies; these latter tXTO features arc not 
particulajiy easy to interpret at this point, but will be discussed 
in the section concerning remnant corrections. I'Jhat should be obvious, 
however, is that the visual gain is considerably smaller (-10 dB) 
than the vestibular gam, over most of the frequency range, excluding 
the very low-frequency crossover region. 


5.3 Non-linear Least-squares Curve Fit Results 


A non-linear curve-fitting program was used to fit the above data 
v’ith different types of specified transfer functions. Shoi<m in figure 
5 is the gain data for with both the mean and standard deviation 

bar-plotted at each frequency. Superimposed on this data are three 
curves resulting from the fitting program and the choice of three 
transfer functions: 
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The complex pole format indicated above was chosen so as to allow the 
fitting program the greatest flexibility in minimizing the fit error. 

As can be seen from the figure, the lag-lead function approximates very 
roughly the gerneral trend of the means, ^7hile the addition of the wash 
out allows for a very good fit at both the Ipw- and mid-frequcncies. 
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Finally, addition of the lead provides for a good fit at the high end 
of the spectrum. Because of the spread in the data (indicated by the 
standard deviation bars), the fit improvement, as one progresses 
from one transfer function to another, is not especially dramatic 
when measured in terms of residual error: the residual is only 

reduced by approximately 10% when the lag-lead is augmented by the 
washout and additional lead. A more dramatic improvement would be 
evident if only the means were used as the data to be fitted, as 
should be clear by an inspection of the figure 

The parameter values obtained from fitting the double second-order 
function of (Ic) are given below: 

Table 1: Farameter Values for CEj|^ Amplitude Fit 


] 

Parameter i 

K 

to 

1 

1 

w 

2 

! 

C ! 
2 


Parameter Value 6 

10.8 

0.29 

0.83 

3.66 

0.85 


Parameter S.D. Og 

1.70 

0.02 

O.U 

0.25 

O.W 



Note that both the poles and zeroes are complex, so that the fit 
procedure does not allow for a simple justification of the functional 
form by appealing to the idea of cascading simpler first-order 
transfer functions This aspect of the problem will be discussed 
later. 

Similar transfer function fits may be made to the "visual" channel 
transfer funct5.oii data specifying these are shown in figure 6. 

The functions used were: 


























Integrator plus lead ; 
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I (x''s + 1 ) (T S + 1 ) = I (S^ + 25 w S + (JO^ ) (le) 

1 4 S C 02 2 2 2 

where the parameter choxce was made for consistency with the parameters 
used in the fits on CE^^. As can be seen from the figure, both functions 
fit the data means quite x^e]l, with the additional high-frequency lead 
allowing for a better fit at the high end of the spectrum. Again, 
because of the data spread andicated by the standard deviation bars, 
the reduction in residual fit error is quite small with the additional 
lead term, but it is clear that the mean trends are better fit with 
the additional Dead. 

The parameter values obtained from fitting the integrator plus 
double lead function of (le) arc given belo\^: 

Table 2. Parameter Values for CE 2 Amplitude Fit 


Parameter 

K 

T 

X ! 

w 1 

5 




1 

4 
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Parameter Value 0 

0.16 

1.26 

0.13 ■ 

2.45 

1.71 


Pai'ameter S.D. cfg 

0.02 

0.04 

0.01 

0 30 

0.30 



Integrator plus 
double lead 


Because of the large frequency difference seen in the break frequencies 
of the data of figure 6, the double zero consists of two real zeroes. 
The complex parameter values (w^, 5 ^) are included here for convenient 
reference later. 

Wiat souDd he clear at this point is that, if the "vestibular” transfer 

fA- 



function is given by (Ic) and the "visual" transfer function 
CE2 by (le) , and the parameter values are as specified in tables 1 
and 2, then a problem arises because of the non-congruence of any 
poles or zeroes. That iSj if the control strategy C is anything 
but a DC gain, then we should expect to see its poles and zeroes 
common to both CE^ and CE2* But this is not the case, which suggests 
that either C is a pure gain so that we actually have measures of 

Ej^ and E2, or that the curve fitting just described is a premature 
exercise, which should await further processing of the data to 
account for control strategy dynamics. This question will be an- 
swered in the next section; for no^7 however, some additional 
observations concerning the phase characteristics are appropriate. 

Shown in figure 7 and 8 are the calculated phase angles associated 
with CE^ and CE2, respectively. As before, the data is summarized 
by bar-plots of the population mean and standard deviataon at each 
frequency. The dashed curve in each figure is the phase predicted 
by the transfer functions just described used to fit the amplitude 
data. Ihe large discrepancy between this predicted phase and the 
actual trends of the data is presumably due, in part, to operator 
latency, as we have seen 30 earler phase data fitting cxercases. 
Addition of a dead-time factor to the transfer functions results 
in the solid curve shown in each figure. The tv?o delay times 
were chosen by visual inspection, and provide an approximate fit 
to the phase data means. 



Two brief cotranents are appropriaio here. First, from figure 7, 

It should be clear that the curve fit to the phase data at low 
frequencies is inadequate in the case of the "vestibular” channel. 
The model predicts a lead^'vdiere we actually observe a lag, and 
would suggest the incorporation of a low-frequency lag term an 
the model. How to accomplish this without affecting the AR 
curve fit IS at present unclear. The second point concerns the 
relatively large delay times necessary to account for the observed 
high frequency lags, approximately 0 6 to 0.7 seconds. Presumably 
this delay time can be apportioned between the operator's estimation 
and control functions; the balance between the two necessDtates 
a closer look at the operator’s control strategy, and is the 
subject of the next section. 





































6.0 VISUAL FIELD VELOCITY NULLING EXPERIMENT 


As noted in, tbe^ previous section, the sub3ect’s control strategy C(s) 
is embedded in the describing functions obtained from both the dual- 
input experiment (series IV) and the earlier velocity-nulling experiments 
(series I, II, III), Thus, the objective of determining the linear es- 
timator functions and E2 has yet to be met. 

This motivated the design of another experimental sen es (VI^ aimed at 
determining the control strategy, so that its effect could he divided 
out of the results already obtained. The experiment xjas designed so as 
to minimize the. estimation task of the subject, but maintain the same 
controller structure used earlier. Specifically, the task chosen xjas 
a standard human operator visual compensatory tracking task 3 n xdiich 
the subject x^as instructed to null visual field velocity via appropriate 
motions of the control x^heel. The disturbance noise injected into the 
projector drive loop x?as identical to that used in the dual- input ex- 
periment (i.e, , the signal d^) and the control x<rheel polarity v?a& chcinged 
to be consistent x^ith the field nulling task, i.e., a right wheel deflec- 
tion resulted in right field mo Lion. The same plant dynamics as before 
(P(s)) x-rere inserted betx^een x^beel deflection and projector motion 

The trainer remained stationary throughout the task, and the subject 
was Informed of this prior to the experiment. In addition, to avoid any 
possibility of circularvection induced by the vd sual field motion, the 
side windoX'JS of the trainer x’ere made opaque and the moving stripe pat- 
tern used earlier Xias projected on the traiixer’ s front xjindoxa. No sensa- 
tions of self-motion x^ere induced by this arrangement, as indicated by 
post-test questioning of the subjects. 

^7 
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This experimental serxes can be directly compared to the earlier ve- 
locity nulling tasks in which the visual field remained fixed. Both 
are single input experiments (visual field veHocity in the present series, 
trainer velocity in the earlier ones) both requiring a single estimator 
cascaded with a control logic. Here, hoxi/ever, because of the kno\ra 
high-pass characteristics of the visual system it is reasonable to 
assume that the visual field velocity estimator necessary for this task 
has essentially unity ga^in over the frequency range of interest (non- 
linear gain characteristics are another matter, hoxaevcr, and \7ill be the 
subject of another study). Thus, what is measured in thDs experiment 
is the subject's control logic C(s). 

. Of course, it can be argued that the subject's control strategy in this 
experiment will differ from that used in the motion nulling experiments, 
simply on Uie basis that the tasks are different. Hox^ever, this can be 
countered by noting that the same control whee3 , plant dynamics , and 
input disturbance (in the visual channel) are used in both experimental 
series. In addition, the tasks are similar in that a nulling of per- 
ceived velocity IS being asked of the subject, in one case self -velocity 
and the other, visual field velocity. Since there is obviously no con- 
clusive way to disect the control strategy from the estimators proposed 
in the parallel channel model, it seems reasonable to assume an identity 
between the control logics of the two tasks, especially in view of task 
similarities . 

Shown in figure 1 is a block diagram of the nulling task, summarizing 
the basic features of the experiment and the conjectured functional 
structure of the human operator. From this diagram and the discussion 
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concerning cross-correlatxon functions given earlier , it is a direct 
matter to show that the control strategy C(s) is defined at the input 
disturbance frequencies ^he following relation* 


KC(m^Uo) 




( 1 ) 


, (in W ) 
^tJ2d2^ 1 0 


W (iD.W ) 
2 1 0 


SO that direct calculations from FFT results may be utilized 


In this experimental series, 6 subjects attempted to null field velocity 
for two full periods of the disturbance signal (T = 128 sec.)> for a total 
individual run time of 256 seconds. The FPT’s were performed on each 
128-second segment, so that t'oo estimates of KC were obtained for each 
subject. It IS appropriate to note that 4 of the subjects participating 
in this series also participated in the DI series (IV) , so that it will 
he possible later to divide out control strategy, on an individual-by- 
individual basis 

Shown in figure 2 are t\?o sets of amplitude spectra obtained from one 
subject, illustrating the relatively large compensatory response at the 
input frequencies (circles) and the small remnant (dots) . As was done 
earlier, a smooth curve has been draxim through the remnant amplitudes, to 
provide a simple approximation for remnant corrections to be introduced 
later. A comparison with similar data obtained from the self-velocity 
nulling experiments (senes I through IV) shox<rs considerably higher 
response amplitudes for the present task, indicative of either a higher 
gain in the control block C, or of attenuation In the self-velocity 
estimator blocks (E^ and E^) . This is also reflected in the higher S/N 
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ratios seen rn this task, In this instance the computed S/N power 
ratios arc 34.8 and 33.8 dB for the sub:) act’s two presentations, (to 
be compared with considerably lower 8/N ratios for the DI experiment, 
for example) . 

Shoxra in figure 3 are amplitude ratio (AR) population averages for 
the SIX subjects, computed for this task from (1), wjth the bars 
indicating one-sigma spread. Superimposed on these data are corres- 
ponding AR's which have been corrected for each individual’s remnant 
power, so as to get a more accurate estimate of the linear transfer 
function C(s) , The remnant correction method will not be discussed 
here, as it is the subject of a later section, t^hat is important to 
note, however, is that the two data sets, corrected and uncorrected, 
are obviously insignificantly different, and thus the simple traxisfer 
function approach of (1) is quite adequate for obtaining an estimate 
of C(s), especially in view of the data spread. 

From the figure, the gam trends are quite obvious- integration at 
the lot/ frequencies, followed by a mid-frequency lead break, similar 
to what was seen in the AR data defining the "visual" estimator/con- 
troller function CE^- This motivates the fitting of the following 
function to the AR data: 


-(Ts 1) 
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( 2 ) 


The dashed line in the figure shows one such fit obtained by a visual 
fit of asymptotes; the parameter values are given by: 


K = 1.8 


T - 0.45 


(3) 
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It should be noted that these are not the result of a least-squares 
fitting algorithm (v7hich will be done shortly) , but simply rough 
estimates. It is also appropriate to note that the lox? frequency 
gain drop-off in the data is greater than the 20 dB/decade attribu- 
table to a simple integrator, and the mid-frequency break seen in 
the data is sharper than a break associated with a simple first-order 
lead term. These discrepancies betXvTeen model and data may motivate 
the use of higher-order functional curve fits; for nox<r, hox^ever, the 
basic integrator plus lead function appears adequate. 

The corresponding phase data obtained from this experimental series 
are presented in figure 4, and although not fully corroborative of the 
gain data, do shox? a constant phase lag of 120° at the low frequencies 
(where x\»e xrould expect a 90° lag xJith a simple integrator). Also ske- 
tched on the f3gure is the predicted phase curve for the transfer func- 
tion defined by (2) and (3) . The discrepancy betx^een the curve and the 

data at high frequencies is presumably due to human operator latencies , 

S 

which can be modelled by a dead- time terra, e d . The resulting phase 
shift is seen in the second curve on the figure, a reasonable approxi- 
mation to the phase data is obtained by choosing a dead-tjme of 0.32 
seconds. Again, a better fit can be expected x^iLh the eventual ap- 
plication of a least-squares program to the data 

As noted in the beginning of this section, the plots defining C(s) 
can be used to infer the estimator functions and E2) from the DI 
experimental data. One approach, given in the next section, is to 
simply divide 6-sub3ect means from the txjo experimental series; that 
is, at each frequency, divide the mean value for jCE^j^j by the mean 














value for |c|, to obtaxn an estxmate for 1e^!> sxmxlarly for 
]E 2 j* Obvxously, the same approach can be used for the phases. 

An alternatxve approach, also presented in the next sectxon, xs to 
dxvxde each xndivxdual's DI data by bxs manual control data, and then 
fxnd average values for the resultant estxmates of and E^. It 
was noted earlxer that only four sub:]ects partxcxpated xn both experx- 
mental serxes, so that the populatxon xs a subset of what has been 
considered so far. To show that the populatxon. average Bode plots 
for the four-subject populatxon are not sxgnxfxcantly dxfferent from 
those obtained for the sxx-subject populatxon, xt is only necessary 
to look at fxgures 5 and 6., whei'e the two data sets are superxmposed 
on one another, A formal test of equivalence involves t-tests on 
curve fxt parameters, but curve-fittang of each data set is yet to 
be done. It should be obvxous however, that the two data secs are 
not signxfxcantDy dxfferent, and thus we mxght expect that the two 
approaches used xn the next sectxon to estimate the esLxmator func- 
tions will yield essentxally the same results 
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7.0 ESTIMATOR FUHCTTONS FOR PARALLEL CHARNEL MODEL 


In section 5 , we obtained the gaxn and phase p3ots defxning the 
combined controller/estimator functions for the parallel channel model, 
CE^ and CE^, and in the previous section we obtained similar plots which 
define the control function C. It is the purpose of this section to 
divide out the control la^7 from the earlier results to obtain estimates 
of the estimator functions E^ and E^. 


7,1 Population Division 

The simplest approach to the problem is to simply divide the six- 
subject means from the two series (IV and VI) to obtain an estimate of 
the population mean for each E^. Thus, at each frequency, divide the 
mean value for [CE^j, obtained from Figure 5 of section 5, by the mean 
value for (c|, from Figure 3 of section 6, to obtain an estimate of the 
mean of • The same can be done for E 2 » and obviously a similar 

approach is applicable to phase angle calculation for both functions. 

The results of such calculations are sho\m in Figures 1 and 2. In 
both figures the closed circles denote E^ values, while the open circles 
denote E^ values, the dashed lines only indicate trends with frequency 
and are not to be mistaken for fitted curves. Several points should be 
noted. First, the gain curve associated with exhibats a rapid drop-off 
at low frequencies and a levelling out at the higher frequencies, behavior 
which is entirely consistent with our notion of a "vestibular" transfer 


























function. Had the experimental equipment allowed for a higher range 
of test frequencies to be used, we presumably would have seen an eventual 
drop-off in the gam, corroborating the previous findings of vestibular 
dynamics characterized by a bandpass filter centered around the frequency 
regime of normal physiological motions. The phase curve associated with 
is shown in Figure 2 and illustrates the expected 90° lead at loxj 
frequencies (characteristic of a "washout") with a gradual drop-off at 
mid-frequencies. With a simple washout we would expect zero lag at the 
highest frequencies; presumably the lag seen in the figure is due to 
operator latencies. 

The behavior of the data characterizing the visual channel transfer 
function E2 sho^T^s a considerable contrast. As noted in an earlier section, 
the visuai gaan (Figui'e 1) is a good deal than the vestibular gam 

over much of the frequency range, T7ith approximate equality occuring only 
at the very low end. In this regime, the visual gain is approximately 
constant with frequency, in agreement with our knowledge of circular- 
vection a non-zero gain at zero frequency. Hon/ever, the gain increase 
seen in the mid-frequency regime is unsettling, since it suggests increased 
velocity sensitivity ■^-'ith higher frequencies, behavior which is in conflict 
^7ith the complementary filtermg model. In fact, \7hat we would have expected 
is a drop 5n gain with increasing fiequcncy, so that the only source of 
high frequency motion information would be by the vestibular pathx’ay. 

The data shov?n here contradicts such a notion. 

There are at least four points worth considering before arriving at 
any definite conclusions regarding the results 311st presented. The first 
Is that the results were obtained by dividing the data from one six-subject 



population by the data from a different six-subject population. Intcr- 
subject variation could be a subtle cause of the unexplained visual 
gam behavior. To resolve this point, transfer function divisions have 
been made on an individual-by- individual basis for the 4-person subset 
of subjects participating in both experimental series. The results are 
the subject of the next subsection. 

The second point is that the resuDts shown here are based on data 
which was not corrected for operator remnant contributions. It maj^ be 
recalled from section 6 that the corrected and uncorrected gains assoc3 ated 
with the control strategy C did not significantly differ, and thus it might 
be surmised that such a correction would be of little value here However, 
it should also be recalled that the visual field velocitj' nulling task 
(Series VI, used to determine C) was characterized by reJatively high 
S/N ratios. This is to be contrasted to the data obtained from the dual- 
input experiment (series IV), m which considerably lower S/N ratios \7er.e 
obtained for the visual channel (recall figure 2 of section 5). Thus, one 
might expect that the visual transfer function inferred from the DI 
experiment, CE^, would shon? a more significant change ^dien remnant 
corrections are incorporated. Presumably, such a change would be reflected 
in the calculated gain and phase for presented here, and thus "correct" 
the increasing gain tendencies seen in Figure 1. The results of such an 
approach will bo discussed in section 8. 

The third point reflects on the validity of the proposed parallel 
channel linear model: the model falls to incorporate a non-linear gain 

characteristic, a property often associated with human operator behavior. 
Specifically, we note that for the visual channel, the gain increases 



at the hxgher frequencies. But from Figure 3 of section 4, we should 
recall that the amplitudes of the input disturbances decrease with 
increasing frequency. Thus, the possiblity exists that what we are 
observing is not a gain dependence on frequency, but on amplitude; a 
gam increase with decreasing amplitudes . This type of non-linear gain 
function IS sketched in figure 3, and is consistent with behavior seen 
in subjective estimation performance in psychophysical experiments 
Consideration of this problem is a current area of effort. 
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A final point ^70rth noting is that the parallel channel model, as 
IS, may he a quite valid representation of the human velocity estimator, 
and the data trends of figures 1 and 2 may be an accurate description of 
the two 'individual transfer functions of the model. What then remains 
to be resolved, however, are the conflicts between the predictions of this 
model and the results obtained in past circularvectiou studies. For instance 
if the visual channel does indeed have a significant gam at high frequencies 
why then does the observed clrcularvection sensation build up so slowly, 
in response to a step input of visual field velocity? Or, more to the 
point, why do high frequency (say 1 Hz) oscillations of the peripheral 
visual field fail to elicit a circu3 arvection sensation, whereas low 
frequency (say 0.01 Hz) oscillations do? 



7.2 Individual Divjsion 


Shown in figures 4 though 7 are gain and phase plots showing average 
values for and calculated by dividing an Individual’s response in 
one experimental session (IV) by his response in the other (VI), and then 
averaging over the 4-subject test population. This approach corrects for 
intersub3ect gam variation, and further, allows for a calculation of gain 
and phase standard deviations at each frequency. One sigma bars are showui 
on the figures. 

A comparison of these results T^xth those described above (figures 1 
and ?) shoii’s that this method of transfer function estimation yields 
essentially the same results, as expected. It may be noted that the 
gains are generally lower here than those of figure 3, and that the lo\? 
frequency phase leads for both functions arc lower than those of figure 2, 
but, otherwise, the results are basical] y unchanged. Specifically, the 
increase in the visual channel gam ^^ith frequency is still present, and 
cannot be accounted for on the basis of intersubject variations. 

Tile figures also illustrate proposed transfer function curves, obtained 
from a simple visual inspection of the data: a washout for E^, and a lead- 

lag for E2* The actual parameter values are only rough estimates, and not 
the results of a least squares curve fit program. It is evident that both 
functions do a fair job of fitting the AR data means; the phase data , 
however, requires the addition of dead-time terms lo follow the phase 
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lag trends at high frequencies. Although this type of latency can be 
justifxed by assuming some type of estimatxon computation time, it should 
be clear from figure 7 that the augmented transfer functxon fit to the 
visual channel phase is not particularly satisfying, especxally in the 
mid-frequency range. 

At thxs poxnt, xt xs appropriate to investigate the effect of remnant 
corrections on the data, to see if the visual channel gain trends are 
affected, and if such corrections afford an improvement in transfer 
function curve fitting. 

Part 8 of this section of the report is m preparation. 





